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Abstract 

Physical fractals invariably have upper and lower limits for their fractal struc- 
ture. Berry has shown that a particle sharply confined to a box has a wave 
function that is fractal both in time and space, with no lower limit. In this ar- 
ticle, two idealizations of this picture are softened and a corresponding lower 
bound for fractality obtained. For a box created by repeated measurements (a 
la the quantum Zeno effect), the lower bound is Ax ~ At(h~/mL) with At the 
interval between measurements and L is the size of the box. For a relativistic 
particle, the lower bound is the Compton wavelength, Ti/mc. The key step in 
deriving both results is to write the propagator as a sum over classical paths. 
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I. INTRODUCTION 



Confinement in a box with hard walls makes things rough. Berry |I| showed that the 
wave function of a particle so restricted is nowhere differentiable, neither in time nor in 
space. Specifically, let a particle be confined to a box, B, in which its dynamics are free 
but for which its wave function, t/>, has Dirichlet boundary conditions on the walls. Let i/j 
initially be constant in B. Then as a function of time the graph of tp (e.g., of Reip) at a 
fixed point has dimension 7/4. For fixed time, the graph of ip is a hypersurface of dimension 
D = d + 1/2 where d is the spatial dimension of B. The boundary area of B is assumed 
finite. (Other well as nongeneric exceptions are also covered in PJ.) 

In a recent article 0, Facchi, Pascazio, Scardicchio and I found a way to build a box 
using the quantum Zeno effect (QZE). You "measure" at intervals At whether a particle 
is in a region, B. As At —>■ we showed that the particle does not escape, which is the 
standard QZE. Our main point was that within the box there is free particle dynamics, with 
Dirichlet boundary conditions on the walls of B. 

Both HI and || idealize unnaturally: you don't find perfect boxes in Nature, you don't 
create wave functions with sharp discontinuities and you don't measure infinitely often. In 
this article I will show how backing away from one of the idealizations softens the extreme 
predictions of fractality. Put differently, observed physical fractals generally have upper and 
lower ranges of validity. What I show here is that when the quantum Zeno box is made 
penetrable by virtue of finite time intervals between measurements, the fractality softens 
below a certain distance scale. Specifically, fractality is lost when 

Ax < At— , (1) 
mL 

where At is the interval between measurements, m is the particle mass and L the box size. 
This relation is consistent with the results of [§,[§], in that a particle confined to a box the 
size of its Compton wavelength travels in straight lines. 

This leads me to examine a second implicit idealization. Both ]I| and |2j use the 
Schrodinger equation as a descriptive framework. The path integral of this equation uses 
Brownian-motion-like paths, rather than those smooth on the Compton wavelength level, 
as is the case relativistically. Interestingly, it is possible to recover this same Compton 
wavelength cutoff of fractality starting from our path integral solution of the Schrodinger 
equation, plus the existence of some maximal velocity. 

The plan of this article is first to recall the results of on how to build a Zeno box, 
and how to let it leak. Then I will rederive Berry's results in a heuristic way using path 
integral language. These perspectives will then be combined to show how the leaks in the 
box truncate the lower scale for fractality. Finally I will show a second kind of fractality 
breakdown, arising from limitations of speed. 

II. THE ZENO BOX 

A particle of mass m is put in a one-dimensional box of length L and at intervals of time 
At it is checked to see whether it is still in the box. The result of is that as At — ► 
the particle remains in the box and within the box evolves under the free Hamiltonian 
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with Dirichlet boundary conditions. This assertion is proved by studying the (combined) 
propagator for time-t free evolution, followed by projection: 



G(x,P,y) = (x\E B U(t)E B 

= Xb( x ) 



m 



2nith 



exp 



X B (x)(x\U(t)\y)x B (y) 
im(x — y) 2 



2M 



x B (y), 



(2) 



where Eg is the projection operator on the box, x B ^ ne characteristic function of the set 
B (so that E B applied to a state is multiplication by X B { x ))i an d U(t) is the operator of 
unitary time evolution, exp(—iHt/h), with H = p 2 /2m. 

The use of the Trotter product formula to analyze Eq. (^]) is inappropriate. Considered 
potential, logx s (^) is singular. 

The method of was to express G in a basis of eigenf unctions of the Dirichlet Hamil- 
tonian (a selection that does not preordain the conclusion) and develop it in an asymp- 
totic short-time expansion. For the box, B, defined by the interval [0, L], the basis (with 
Hu n = E n u n ) is 



U n [X) 



2 ( nirx\ 
sin I I 



L~~\ L 
The matrix elements of G are 

G mn {t) 



E n 



h 2 n 2 n 2 
2mL 2 
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1,2, 



dx l dy u m (x) ., 
o Jo y ' V 2-nith 



m 



exp 



im(x — y)' 
2kt 



u n {y) ■ 



(3) 



(4) 



Small-t asymptotic analysis of (ffl) yields 



it 



G mn (t) = 5 mn (l--E n )+0(t 3 / 2 ) 



h 



(5) 



For t — > this recovers the free propagator with Dirichlet boundary conditions. For the N- 
projection QZE formulation, showing that the remainder in Eq. (j5|) is 0(t 3 / 2 ) is sufficient to 
show that for t = 1/N, and N projections, one obtains lossless confinement. Moreover, the 
fact that this turns out to be diagonal in the Dirichlet basis shows this to be the self-adjoint 
extension of H chosen by the Zeno projection process. Eq. (]5|) can also be used to study 
leakage, that is, imperfect confinement for finite measurement frequency. 

In [0 this result is extended to many dimensions, to complicated boxes, and to Hamil- 
tonians with potentials. But these matters do not concern us here. 

The point that I do take is that to create a hard wall, one does not need to be there all 
the time, just very frequently. 



III. FRACTALS IN A BOX 



Before getting fractals, we will see how not to get fractals. Although the initial wave 
function ip(x, 0) ee ■^wX[o,l]( x ) is discontinuous (constant in the box, zero outside), if it is 
allowed to propagate freely, that is, on the entire line, it is perfectly well behaved. 

Using the free propagator (on the entire line), at time-t this wave function becomes 
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V>OM) = -7= / dyg(x -y,r) , where g{u,s) = - 7 == and r = — . (6) 



This is clearly a smooth function of x. We also find it convenient to look at the derivative 
of if), in particular 



dtp 1 



Ox yfjj 



g(x,r)-g(L-x,r) , (7) 



a result derived by changing d x to d y within the integral. Again, this is smooth. 

Returning to fractality, Berry's arguments depend on the following property of slowly 
converging Fourier series. Consider 

f(u) = V a m e mu , with \a m \ 2 ~ -. — for \m\ — > oo and 1 < (3 < 3 . (8) 

— ' 777 " 

in I ""I 

If the phases of {a m } are random or pseudorandom, then it follows that the graph of Re/, 
Im/, and generically \f\ 2 are fractals of dimension 

D=±(5-(3). (9) 

I will give a heuristic proof of this result, based on a criterion immediately derivable from 
box-counting definitions of fractal dimension ||. If the graph of f(u) has dimension D, its 
mean square increment over a small distance Au satisfies 

((A/) 2 ) ~ (Au)"' 20 . (10) 

The square deviation of / is given by 

'mAu\ „ . / nAu" 



\f(u + Au)-f(u) 



£ £ a m ay u ^2 sin (^) 2 sin 



(11) 



The average is performed by invoking the "pseudorandom" phases to eliminate all n ^ m 
terms in Eq. ([11]). This leads to 



«A/)»> = 4Ek|W „ ^Irf (=£) (12) 

with K a generic constant. For small Am, and up to small errors and irrelevant factors this 
can be replaced by an integral and rescaled: 

Combining Eqs. flT3"|) and (fl(]| ) one obtains the desired result, Eq. ([|). 
The propagator within the box is given as usual || by 

_ 2 /rnrx\ f mxy\ ( hn 2 Tr 2 \ 
G(x, t; y) = E Z sm (-^-J sm j exp ^ J (14) 
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Recall the definition of the Jacobi theta function |7|], 



6 3 (z,s)= ^2 ex P (insn 2 + 2iznj . (15) 



Shuffling terms in Eq. ([Tj) yields 



^i/ s 1 „ ( 7T . . 7T^t \ 1 . ( 7T - . irht \ . . 

G <^ (; ») = 21* (si* 1 - »>■ "SO?) - 2Z* + »)• J (16) 

Several properties of #3 now play a role. First it has period ir and quasiperiod its, 

e 3 (z + 7r,s) =9 3 (z,s) , and 3 (z + 7rs,s) = e"^" 2 "^^, s) • (17) 

There is a zero of # 3 at z = (1 + s)tt/2, hence, by Eq. (|17|) , there is an infinity of them. 6 3 
also satisfies an identity of far-reaching significance 

9 3 (z, s) = (-isyWef^'Os -1) . (18) 

If Im s > 0, #3 is clearly an analytic function of z. On the line Im s = things get interesting 
and in fact, from Eq. (0), this is the case for the physical propagator. In particular, the 
countable set of zeros now all lie densely on the real line — this for an integral kernel that 
preserves norm. I have been fascinated by this peculiarity of the propagator for a long time 
||, and am not surprised that the wave functions propagated by G can be rough. 

Applying the fundamental identity, Eq. (|18|), to G is straightforward, and a little algebra 
yields 

G(x, t; y) = Y<{g(x -V- 2nL, r) - g{x + y - 2nL, r)} (19) 

n 

with the same g and r that were defined in Eq. (^). This shifts the propagator from an 
energy expansion to a path expansion, since each term in Eq. (|l~9| ) can be identified with 
a classical path. The direct path from y to x, not bouncing off any wall, is the positive 
n = term in (|19D. Consider a path from y to x that bounces off the right wall then the 
left, then the right again, and then reaches x. The distance it covers is: L — y (getting to 
the right wall) plus 2L (right to left to right again) plus L — x (getting from the right wall 
to x). This gives 4L — x — y, the negative n = 2 term above. This can be visualized by 
making copies of the box along the real line, with each bouncing path corresponding to an 
origin point at a displaced image of y (rather like the universal covering space picture in 
0). The propagator, (|T9|), could also have been derived by the method of images. Thus 
each term in flT9D satisfies the Schrodinger equation and the particular combination taken is 
that guaranteed to satisfy Dirichlet boundary conditions. 

We next establish Berry's result using the path expansion of the propagator. Again we 
consider a wave function that is initially 1/vL on [0, L], zero elsewhere. Now, however, its 
motion is confined to the box and it evolves under the propagator (|19D . In evaluating the 
spatial derivative of the wave function at time-t, the identity that took us from Eq. (ffi) to 
Eq. ([7D is still valid (although for the "x + y" contribution the sign is opposite), and we get 
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dip 



7 f = —Y^ {g(x - 2nL, r) - g(x - (2n + l)L)} = "£(-l) n g(x - nL, r) 

OX a/ L n n 



(20) 



At this point the reader may be feeling uncomfortable at the display of the derivative of a 
function known not to possess one, but we sidestep this problem by adding a small negative 
imaginary part to any convenient quantity, for example r. Everything is now analytic and 
our formal manipulations justified. 

Our aim is to calculate (|A?/;| 2 ) and use Eq. fllCf ). We write 



Aip 



- 7 ='£(- i ri>^-nL,T) 

V L n 



nL, t)Ax 

g(x + Ax — nL, r) — g(x — nL, r) 
dg(x — nL, r)/dx 



(21) 



Recalling the definition of g, this becomes 



Aip 



2t 

Wl 



^ ' x-nL 



[x 



exp i- 



nL)Ax .(Ax) 2 " 



2r 



(22) 



The next step is to take the absolute value squared of Aip, written as a double sum, say 
over n and m. We then argue that "random phases" make the n ^ m terms drop out and 
find 



sin 



(x — nL)Ax 
27 



(23) 



This is essentially the same sum as in Eq. (0), now with (3 = 2. It gives a fractal dimension 
3/2, which is the correct value. 

Note that in Eq. (p3|) we could let r lose its imaginary part before considering the limit 
Ax — > 0. This is because the series is convergent without r's help. 

What about the random phase assumption? For sufficiently small Ax the contributions 
to the sum come from large n. For the square, there are two sums (over "n" and "m"). 
The terms g(x — niL, r)g(x — nL, r)* will have phases on the scale of (m 2 — n 2 )L/r, whose 
remainder modulo 2tc provides the needed pseudo-randomness. 



IV. LIMIT OF FRACTALITY IN THE ZENO BOX 

With the path-based perspective just developed for box fractality, we can see how leaks 
in the box affect the roughness of the wave function. 

In the Zeno box, the particle is contained through repeated localization measurements, 
implemented at intervals At. The propagator for time evolution to time t (with iV = t/ At 
projections) is the N^ 1 power of the propagator of Eq. (FJ). As shown in ||, in the limit 
N — > oo this goes over to the box propagator, given for example in the path expansion form, 
Eq. ©. 

Now consider the contributions to the sum in Eq. (|l^) or (^3|) . Those with summation 



index n bounce off the walls about n times, so that if the final time at which the wave 
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function is evaluated is t, they are hitting the walls at intervals t/n. But these terms should 
not appear in a finite- N Zeno-box propagator if n > N\ For in that case, there would be 
nothing to stop the particle from passing the often-nonexistent boundary. It follows that 
the only paths that could appear in an approximation for the Zeno-box propagator must 
satisfy n ^. t/At. 

But now return to Eq. (|T3|), a special case of which is encountered in Eq. (^). The bulk 
of this integral comes from contributions for which y = 0(1). Thus, for example, for (3 = 2 
(our case), the range 1/4 < y < 4 includes more than 75% of the integral. Now "y" in 
Eq. ( P^ ) is nAu, and Au in Eq. ( |2"BD is LAx/r. It follows that the fractality condition is 
nLAx/r ^ 1 or r/LAx ^ n. 

Recalling the definition of r (Eq. @), from the last two paragraphs we have the joint 
requirements 

fit t . nj . 

& n and n % — — . (24) 



mLAx " At 
It follows that ip can exhibit fractal structure only for 

Ax > -ArAt , (25) 
mL 

which is the same as Eq. ([!]), but stated as a converse. 

There is a second way to arrive at this conclusion. Instead of the fractal dimension for- 
mula, we make use of the fundamental identity, Eq. (|18|). As above, the sum-over-classical- 
paths expansion, ([B]), leads to the condition n ^ t/At for the Zeno box, with At the interval 
between measurements/projections. But we now implement this restriction in the propaga- 
tor (|19D by means of a cutoff — a gradual one, achieved by making t complex |J. Recall that 
each term in Eq. ([HJ) is of the form [K/^/i] exp (im(£ — 2nL) 2 /2ht), with K a constant and 
£ = x ± y. Letting t — > t — it with e > introduces the real factor 

f-em(£ - 2nLf\ f-2emn 2 L 2 \ , , 

CXP ( 2W ) ~ CXP ) ■ < 26 » 

To cut off the sum in (^) at some iV (= t/At), set the argument of the last exponent in 
Eq. ( PB[ ) to 1, leading to 

M 2 

(27) 



2mN 2 L 2 



The next step is to go back to the eigenfunction expansion by means of a second application 
of Eq. flTS|). The result is just the formula ([14]), except that t is complex. In particular the 
argument of the exponent giving the time dependence acquires a real part: 



. Tin 2 it 2 ^ 
2mL 2 



exp —it 



( hn 2 n 2 \ ( n 2 TT 2 h 2 t 2 \ 

exp \-^^l 2 ) =exp r^wj • (28) 



Note the distinction: the lower case n is now the index in the eigenfunction expansion, while 
iV is the number of measurement projections. From Eq. fl28|) we see that the cutoff in the 
path summation around N leads to a cutoff in the eigenfunction expansion around 
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2mL »r 

n = — - — N . 29 

Eigenf unctions with greater n will have reduced amplitude. Recalling the form of the eigen- 
function, sin (n7rx/L), we see that the cutoff will suppress detailed variation of the wave 
function on a scale Ax ~ L/mr. Substituting the value of n given in Eq. (|29|) yields 

fit h . 

Ax — = -At . 30 

2mLN 2mL v ; 

Aside from an irrelevant factor of 2, this is the same as our earlier result. 



V. RELATIVISTIC LIMITATIONS 

For relativistic particles it is known that the agitated Brownian-motion-like particle 
paths characteristic of the path integral smooth out into straight lines. The path expansion, 
Eq. (|19"D, will allow us to see how this destroys fractality. 



A. At the Compton wavelength scale 

A particle of mass m has a characteristic length, its Compton wavelength, A c = Ti/mc. 



Rewriting Eq. (f25[) in terms of this quantity yields 

Ax ^ At 



> 



A c L/c 



(31) 



for the minimal level of fractal structure. For the right hand side of this inequality to be 
small, the interval between measurements would have to be small compared to the time it 
takes light to cross the box itself. Even without claiming to understand the measurement 
process, this would seem a limiting time interval, which in turn implies that the particle 
does not have fractal structure below the scale of its Compton wavelength. 

It would thus appear that the QZE plus a causality-like condition are enough to fix A c 
as a lower level of fractality. As we will see, the result is more general and for box-fractality 
is more the result of the path expansion than of the QZE. 

B. The checkerboard path integral 

These considerations bring to mind properties of relativistic particles that are known from 
another path integral, the Feynman "checkerboard" sum for the Dirac equation. In his early 
work [ 1(| , Feynman developed a sum over paths forumulation for the 1 + 1 dimensional Dirac 



equation. This path integral, independently discovered by Riazanov was published by 



Feynman in ]E^. In parallel, a sum over paths for the telegrapher equation was proposed 



by Kac ||13|| . It was later found 0] that this is an analytic continuation of the checkerboard 
sum and that the underlying wave equations (Dirac and telegrapher) bear the same relation 
14| . These formulations have been developed by many researchers over the years [|3"|,f4l,|15 
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have been extended in various ways to 3-dimensions [16], and applied to other problems as 



well |7|. 



In using the checkerboard formulation, we found that at the relativistic level particle 
trajectories become smooth. Instead of the nonrelativistic path integral's Brownian-motion 
paths, fractal objects with Aa; ~ V~At, the particle travels in straight lines at the velocity 
of light, but at random times reverse direction. These reversals are Poisson distributed 
with rate mc 2 /h. This implies that the particle reverses after traveling on the order of one 
Compton wavelength |TJ| . 



The nonrelativistic limit of this, and the emergence of jerky, fractal paths, now becomes a 
standard piece of mathematics — modulo the usual "1" For fixed distances, you allow the time 
scale to grow. Consider the physical Brownian motion of say, pollen. A reasonable model is 
to assume that between collisions the path is a straight line and that the collisions themselves 
are Poisson distributed. Brownian motion then emerges as a limit (and then an idealization) 
of this phenomenon. A quantitative version of this transition can be implemented for the 
nonrelativistic limit of the checkerboard path integral. Our picture, described above, is 
"straight line motion at velocity c, with reversals on the average once per h/mc 2 " — leading 
to an expected distance covered between reversals of h/mc. Now the diffusion constant 
derived from microscopic models of Brownian motion is on the order of D ~ (Ax) 2 /At. 
Putting this together for the checkerboard process, one gets 

B~^ = ^. (32) 
h/mc z m 



From many considerations [15], this is the right diffusion constant to assign to Schrodinger 
particles. 

From 0, the condition for the validity of this procedure is mc 2 /h 3> cd x . This inequality 
comes from a proposed expansion of the propagator, so that the characteristic momentum 
would be m5x/St (with 5x and St increments of the arguments of the propagator). Replacing 
d x by mSx/hSt, the quoted condition becomes Sx <C cSt. 



C. Practality breakdown from the existence of a limiting velocity 



Starting from the Schrodinger equation the same conclusions about smoothing can be 
drawn from the mere existence of a limiting velocity, c. Explicitly the expansion of the box 
propagator, Eq. (|T5j), is 



G(x,t;y) = 



m 



lixith 



exp 



im(x — y — 2nL,y 
2kt 



exp 



im(x + y — 2nL)" 
2ht 



(33) 



For large n, the rfi 1 path has covered a distance 2nL, and has done it in time t. The 
existence of a limiting velocity implies that such terms should not appear if 2nL > ct. In 
other words, we only allow paths with n < ct/(2L) By the same arguments that followed 
Eq. ( pi]) a cutoff in n in the integral ([13]) implies n > htj (mLAx). Combining this with the 
constraint just derived from the existence of a limiting velocity, we have 



ct 
2L 



> n> 



ht 



LmAx ' 



(34) 
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implying 



Ax > — . (35) 
mc 

A limiting velocity implies a limit on fractal structure. If c is identified as physical light 
velocity, the limit is the Compton wavelength. 
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